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We report the experimental measurement of the equation of state of a two-dimensional Fermi 
gas with attractive s-wave interactions throughout the crossover from a weakly coupled Fermi gas 
to a Bose gas of tightly bound dimers as the interaction strength is varied. We demonstrate that 
interactions lead to a renormalization of the density of the Fermi gas by several orders of magnitude. 
We compare our data near the ground state and at finite temperature to predictions for both fermions 
and bosons from Quantum Monte Carlo simulations and Luttinger-Ward theory. Our results serve 
as input for investigations of close-to-equilibrium dynamics and transport in the two-dimensional 
system. 
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The rich phenomenology of fermionic many-body sys¬ 
tems reveals itself on very different scales of energy, rang¬ 
ing from solid state materials and ultracold quantum 
gases to heavy-ion collisions and neutron stars. Under¬ 
standing the underlying mechanisms promises substan¬ 
tial advances both on a fundamental and technological 
level. Ultracold quantum gases provide a platform for 
the exploration of the macroscopic phases and thermody¬ 
namic properties of fermionic many-body Hamiltonians 
in a highly controlled manner [T]. In particular, using 
strongly anisotropic traps, it is possible to enter the 2D 
regime [MZ] which is of large interest to the condensed 
matter community [iin]. 

The thermodynamic properties of a many-body sys¬ 
tem are encapsulated in its equation of state (EOS) 
n(fj,,T, {gi}), which expresses the density n as a func¬ 
tion of chemical potential /i, temperature T, and further 
system parameters {gi} characterizing, for instance, the 
interactions between particles. For ultracold atoms with 
short-range attraction, the only additional parameter is 
the s-wave scattering length a. This universality allows 
one to describe different atomic species by the same EOS 
n{p.,T,a). The equilibrium EOS is also the basis for 
studying dynamics close to thermal equilibrium. 

In this Letter, we report the experimental determi¬ 
nation of the EOS of two-component fermions with at¬ 
tractive short-range interactions in the 2D BEC-BCS 
crossover regime. We tune the interaction strength using 
a Feshbach resonance to connect the well-known limits of 
a weakly attractive Fermi gas and a Bose gas of tightly 
bound dimers. We report the measurement of the finite 
temperature EOS in the intermediate, strongly correlated 
region and compare with theoretical predictions. 

Our experimental setup consists of a population- 
balanced mixture of A^ ~ 100,000 ®Li-atoms in the lowest 
two hyperhne states, which we denote by |1) and |2). The 
interactions between both species can be tuned by means 
of a magnetic Feshbach resonance [Bill]. The atoms are 


trapped in a highly anisotropic trapping potential, which 
is radially symmetric to a high degree in the ccy-plane and 
provides a tight confinement along the z-direction with 
the aspect ratio of frequencies Wa, : Wy : a;2 « 1 : 1 : 310. 
A detailed description of the experiment is given in [7]. 
This strong anisotropy induces a quantum confinement of 
the many-body system, with discrete excitation spectrum 
in the z direction and an effective 2D continuum of states 
in the xy plane. A 2D system is realized if the interacting 
system is in its ground state in the z-direction. Although 
the explicit form of this ground state is not known and 
depends on the 3D scattering length csd [12], the ratio 
N{ujx/ujz)'^ ~ 1 close to unity indicates that the first ex¬ 
cited state in z direction will be partially populated. In 
this quasi-2D regime, interactions can still be described 
by an effective 2D scattering length 020 Pda HD. 

In order to parametrize the local density of the gas, we 
introduce the Fermi wave vector kp = (27rn)^/^, where 
n = ni -I- 772 is the planar density and rio- = n/2 is 
the density of atoms in state |cr). We further dehne the 
Fermi energy ep = /i^fc|/(2M) and Fermi temperature 
Tp = ep/A:B, where h is Planck’s constant, kp Boltz¬ 
mann’s constant, and M the atomic mass of ®Li. The 
2D crossover [T5H24] is parametrized by ln(fcpa 2 D)- In 
the BEG limit ln(A:p 02 D) "C —1 the system can be de¬ 
scribed as a 2D Bose gas with effective coupling constant 
g « —27r/ln(A:pa 2 D)- In the BCS limit ln(fcpa 2 D) 1 
the thermodynamic properties approach that of an ideal 
Fermi gas. In analogy to the 3D BEG-BCS crossover 
[25H^ , the 2D gas undergoes a finite temperature phase 
transition to a superfluid state for all values of the param¬ 
eter (fcpOsD)”^ |7]. The associated transition, however, 
is of BKT type and the superfluid phase exhibits quasi 
long-range order [2M34] . 

For a 2D ultracold quantum gas with short-range at¬ 
traction, a two-body bound state with binding energy 
£b = ^^/(Afa 2 p)) > 0 exists for all values of the cou¬ 
pling strength. This has to be contrasted to the 3D case. 
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where a bound state only exists on the Bose side of the 
resonance. In both cases, the chemical potential for a 
single fermion becomes negative in the Bose limit and 
approaches /i « —ee/S. The crossover point between the 
Bose and Fermi sides can be defined by the zero crossing 
of /i [3S]. The chemical potential shifted by the bound- 
state energy, jx = /r -|- £ 3 / 2 , is positive at high phase- 
space densities. At zero temperature, n > 0 is equivalent 
to /2 > 0. 

In the experimental realization of the quasi-2D gas 
in an anisotropic 3D trap, the interaction strength 020 
depends on the typical momenta of scattering particles 
and thereby on the filling in the trap with axial fre¬ 
quency (jJz- One can write 020 = 
where is the scattering length in the dilute limit 
and Arc is a positive function which reduces 02 D at fi¬ 
nite density O [131 EH] • The correction to vanishes 
in the Bose limit where /i —>■ 0, and becomes strongest 
in the Fermi limit where /i ~ ep. In our experiment we 
have huJz/kB = 265nK, which has to be compared with 
typical values /io = (40 ... 200)nK and T = (60 ... 25)nK 
when going from the Bose to the Fermi limit. As most 
particles are in the center of the cloud, we approximate 
a 2 D and ep by their central values using Aw{flQ/Huiz), 
giving Aw « 0.2,0.9, and 1.4 in the Bose, crossover, and 
Fermi regimes, respectively. 

We extract the EOS of the homogeneous gas from 
the trapped system by using the local density approx¬ 
imation (LDA) which assigns a local chemical potential 
— V(r) to each point r in the trapping poten¬ 
tial V{f) [37]. Since V{f) is known to a high precision, 
the homogeneous density n{fi, T) can be deduced from 
the measured local in-situ density of the inhomogeneous 
system n(f) = n(/ro — V{r),T) once no and T have been 
determined [35] . The extraction of the homogeneous EOS 
from the trapped gas has been applied to both bosonic 
and fermionic systems and successfully compared with 
theoretical calculations |33H3H] . 

Low-temperature EOS .—In order to determine the low- 
temperature equation of state n{pL,T —>■ 0 , 020 ) we ex¬ 
tract piQ from a Thomas-Fermi (TF) fit of the central 
region of the cloud. The TF model assumes locally 
ep = c ■ p, for the central density region. This scaling 
is valid for large phase space densities (PSD) nX'^, where 
At = {2TTh?/(Mk-BT))^^^ is the thermal wavelength of 
atoms. We find that the prefactor c only weakly depends 
on the temperature and fitting range at sufficiently low 
temperature and high densities, which confirms the va¬ 
lidity of the linear relation ep (x p (see the supplemental 
material (SM) [36] for details). We fit c in the intervals 
IA = [0.4, O.Sjupeak and Ib = [0.5, Ijripeak for peak den¬ 
sity Upeak) and define pQ = {Po,a+ Po.b)/‘2' as the average 
value of both outcomes. 

In Fig. we show the low-temperature EOS across 
the 2D BEC-BCS crossover in terms of p/ep = 1/c vs. 



ln(kpa2o) 


FIG. 1: Low-temperature EOS across the 2D BEC-BCS 
crossover. The experimental results are obtained from mea¬ 
surements of the quasi-2D gas at the lowest attainable tem¬ 
peratures, which corresponds to TjTp ~ 0.05 and 0.1 on the 
Bose and Fermi sides. The data points shown as diamonds 
(circles) correspond to measurements in the superfluid (nor¬ 
mal) phase. The solid red line on the Bose side corresponds to 
the mean field formula P/ef = —ri~^j4 with rj = ln(fcFa 2 D), 
whereas the dashed and solid green lines on the Fermi side 
display the non-selfconsistent and selfconsistent Hartree-Fock 
predictions 1/(1 -|- ri~^) and 1 — ri~^ for weakly attractive 
fermions. The orange line is the prediction for the ground 
state EOS from recent QMC calculations m- 

ln(fcFa 2 D), where kp corresponds to the peak density. 
The corresponding temperatures (60-25nK from left to 
right) are very low compared to Tp (1500-300nK). For 
the plot we averaged the TF slope c over 30 images at 
the lowest temperatures for each value of the interaction 
strength. We find c to be weakly dependent on temper¬ 
ature as we increase T/Tp by (40 — 100)%, which is a 
necessary condition for the applicability of the linear fit 
of the central region. The statistical error of p/ep is 10% 
within the whole crossover. We estimate the error due 
to systematic uncertainties resulting from the absorption 
imaging, atoms in noncentral pancakes |62j . magnifica¬ 
tion, and the determination of the binding energy to be 
15% and 13% on the Bose and Fermi sides, respectively 

[ZllMl. 

Our measured low-temperature equation of state con¬ 
nects both perturbative limits of the crossover. In 
the Bose limit we compare our results with predictions 
for bosonic dimers of mass 2M, dimer density = 
n/2, chemical potential p,d = 2p, and thermal wave¬ 
length Ad = At/v^. The interactions between dimers 
can be modelled by an effective 2D coupling strength 
g = with = y/h/{Mu)z) = 0.551/rm 

the oscillator length of axial confinement |3H]. We find 
p/ep = 0.024(2), 0.046(4), 0.12(1) for the data points 
with ln(fcFa 2 D) < —0.71 corresponding to effective bo¬ 
son coupling strengths g = 0.60,1.07,2.75, respectively. 
This is in excellent agreement with the perturbative Bose 
gas formula p/ep = g/(87r) = 0.024,0.043,0.11. Further- 
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more, we verify g « —27r/ln(fcF 02 D) for very small g, 
which is a result of 020 — on the Bose side, where 
the filling correction Aw is small. 

Far on the Fermi side, our low-temperature data is con¬ 
sistently below the Hartree-Fock (FIF) prediction /i/ep ~ 
/i/ep = 1 — (ln(A:pa 2 D))~^- However, the additional error 
due to systematic uncertainties [36] is of the same size 
as the statistical one displayed in Fig. [T] such that we 
find consistency with HF theory within the errors of our 
measurements. Note that an extension of the BCS mean- 
field theory toward the crossover, which works reasonably 
well in 3D [49l [50], would give fl/sp = 1 for all inter¬ 
action strengths in the 2D case and thus clearly misses 
the crossover physics m- The ground-state equation of 
state has been investigated theoretically in 
In Fig. [^we compare our measurements to recent QMC 
simulations of the ground state SUES]- For a compar¬ 
ison of different theoretical predictions of fl/ep at zero 
temperature see, for instance. Ref. [Sl|. Our data at 
small but finite temperature lie consistently below the 
zero-temperature prediction. 

Finite-temperature EOS. —While the temperature is 
constant within each atom cloud, it varies for every in¬ 
dividual realization of the gas (“shot”). In our analy¬ 
sis we therefore determine T and /io from each density 
profile and construct the dimensionless PSD fi{x,y) = 
/(/3/x,/3eB) = riAy and normalized density hi{x,y) = 
Pe-q) = n/riQ for every shot i. Here, no{g,T) = 
ln(l -I- is the EOS of an ideal Fermi gas. Fi¬ 
nally, we average fi and hi over 30-150 shots to obtain 
the EOS with very small statistical error, even though 
the thermodynamic parameters vary from shot to shot. 
The values of T and can be found by different meth¬ 
ods: whereas the TF fit determines /io from the dense 
central region of the cloud, fitting a reference EOS to the 
outer low-density regions gives both T and /io- 

We first summarize the reference EOSs used in this 
work. In the perturbative Bose limit of small g, the 
outer wings are described by the HF formula ndA^ = 
— ln(l — The Boltzmann limit for a gas 

of dimers or atoms, respectively, reads no = = 

and The latter two formulas 

are elegantly connected by the second order virial expan¬ 
sion = ln(l -I- -I- 2626 ^^^ [3SJ jSSJ jSSj. In the 

weakly interacting Fermi limit 62 0, and n^X^ ap¬ 

proaches the EOS of an ideal Fermi gas. On the Bose 
side, instead, Eb becomes large and the fermion fugacity 
g/3it _ ^Piii-e^/ 2 ) jg extremely small, suppressing the first 
term of the EOS. However, 62 = up to exponentially 
small corrections and we recover the bosonic Boltzmann 
EOS Ua- = Hence, the second 

order virial expansion has the correct limiting behavior 
and provides a well-defined reference EOS throughout 
the crossover. 

We apply the HF formula for the perturbative Bose gas 



FIG. 2: Phase space density on the Bose side of the crossover. 
The experimental data points are shown as filled shapes. We 
compare to bosonic theory with effective conpling strengths 
g = 0.60,1.07, 2.75. Open shapes represent the EOS extracted 
from the QMC simulation of the quasi-2D Bose gas trapped in 
an external potential with similar parameters as employed in 
experiment. The dashed curves show the classical MC predic¬ 
tion for the weakly coupled homogeneous 2D Bose gas from 
Ref. |56] extrapolated to large values of g. For moderate den¬ 
sities we find good agreement between all three approaches. 
For large densities, however, experiment and trapped QMC 
deviate from the classical homogeneous result which scales like 
the mean field prediction ridA^ = ^/3/‘d-l-ln(^ndAd —2/?/td). 
This may be due to quantum effects appearing at large g, as 
well as the axial confinement. 

to determine T and /to only for B — 692G where g = 0.60. 
For the remaining magnetic fields, i?[G] = 732 — 922, the 
central chemical potential /xq is determined from the TF 
fit of the central region. The temperature is estimated 
by r = (TV -b Tb)/2, where Ty and Tb are obtained 
from second order virial and Boltzmann hts to the outer 
region, respectively. This choice is motivated by the ob¬ 
servation that Ty and Tb give upper and lower bounds on 
the true temperature for the interaction strengths consid¬ 
ered here. We quantitatively compare different methods 
to obtain T and /xq in the SM [36] . In particular, we 
show that /io obtained from the virial and Boltzmann 
fits agrees well with the one from the TF fit, which also 
supports the validity of the TF assumption. 

Bose EOS. — Figure]^ shows the EOS on the Bose side 
of the crossover. Due to the exponentially large binding 
energy in the Bose limit, the logarithmic dependence of 
the EOS on /Heb in nA|- = f{/3g, /3£b) can be replaced by 
the g dependence in nX^ = F{l3jl,g), where F{x,g) can 
directly be compared with bosonic theory for coupling g. 
The experimental data correspond to g = 0.60,1.07,2.75. 
The plotted curves represent the average F{x,g) from 
approximately 30 individual shots for each interaction 
strength. 

We compare our results on the Bose side to classical 
Monte Garlo (MG) and Quantum Monte Carlo (QMC) 
simulations of bosons in order to understand whether 
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FIG. 3: EOS in the crossover regime shown as the density normalized by the ideal Fermi gas no{fi, T) = ln(l + The 
experimental data points (filled shapes) are compared with the second order virial expansion at low values of /3p (coloured 
dashed lines). The displayed errors are purely statistical, with systematic uncertainties estimated at 13-15% and given for each 
data set in Table IV in the SM [35]. We compare our results with theoretical predictions in the 2D BEC-BCS crossover from 
Luttinger-Ward theory (|23|, solid black lines) and fermionic QMC simulations (|2^, dotted black lines), with the corresponding 
values of Peb labelling each curve. Note that the vertical scale differs by a factor of 10 in each panel. 


the fermionic system can be described purely in terms 
of bosons. The classical MC computations are valid for 
the weakly coupled, homogeneous 2D Bose gas in the 
fluctuating regime [SSI [57]. In our case g is large, and 
quantum effects are expected to modify the result. Even 
so, the critical temperature for ^ < 3 is well described by 
extrapolating MC to large g [Ml EHl |59| . We also ana¬ 
lyze density profiles obtained from QMC simulations of 
a trapped Bose gas with similar trapping parameters as 
in the experiment; see Refs. miMiiii] for details. We 
extract the temperature and chemical potential of the 
QMC density profiles from the low-density Boltzmann 
regime and apply LDA to obtain the EOS F{(3jl,g). 

We find good agreement of our measurements with the 
EOS extracted from the QMC profiles. The deviations 
between both EOSs are within the systematic errors in 
the determination of n, T and go and we conclude that 
the fermionic system with g < 2.75 is well-described by 
a strongly coupled quantum gas of dimers. Both experi¬ 
ment and QMC are, however, well below the classical MC 
predictions for large jPg and mean field theory. There are 
two effects which could explain this behavior: On the one 
hand, quantum fluctuations become important for large 
g and high densities. On the other hand, both experi¬ 
ment and QMC are performed in a quasi-2D setting with 
nonzero extent in the z-direction. 

Crossover EOS .—Figure shows the EOS in the 
strongly correlated crossover regime between the Bose 
and Fermi limits. The EOS h{x,y) = n/riQ is sampled 
over approximately 150 shots for each of the magnetic 
fields B[G] = 812,832,852,892. We compare our re¬ 
sults with theoretical predictions for the homogeneous 2D 
BEC-BCS crossover from Luttinger-Ward (LW) theory 
[23] and fermionic QMC simulations [Ml- Our compar¬ 
ison covers a substantial renormalization by two orders 
of magnitude in the density nJnQ. It reveals a maximum 
in n/riQ characteristic of the density driven crossover in 


2D [23]. We find that the maximum of height is 

reached at jPg ~ —/3eB/2-|-ln(2) for large /Seb- The origin 
of this scaling can be understood from the virial expan¬ 
sion in the Bose limit: « 2 exp(2/3/x) = 2 dX jPg — 0, 

which implies n/no « 2/ln(l -|- « 26^^^®/^ at 

g = — £b/2; see the SM [36] for details. The difference be¬ 
tween the LW and QMC EOSs lies within our systematic 
errors from the T and go determination and thus cannot 
be resolved with the present analysis. In the Fermi limit 
where go/fujJz ^ 0.75 is largest, we observe that the fill¬ 
ing correction Aw shifts the EOS slightly upward. This 
effect is minimized for small particle numbers. In a re¬ 
cent work by Fenech et al. [45] the EOS in the normal 
phase for jPe-Q < 0.5 has been determined using ®Li atoms 
in the 2D regime. 

In this work we have measured the EOS of ultra¬ 
cold fermions in the BEC-BCS crossover in a strongly 
anisotropic confinement. Our results connect the pertur¬ 
bative Bose gas, the strongly interacting Bose gas, the 
strongly interacting fermionic superfluid in the crossover 
regime, and the perturbative Fermi liquid as we tune 
the effective 2D scattering length a 2 D using a Feshbach 
resonance. Our EOS data covers both the low and in¬ 
termediate temperature thermodynamics of the system. 
We compare with bosonic and fermionic quantum many- 
body theory and find a remarkably strong renormaliza¬ 
tion of the density n/no- These results provide a basis for 
phenomenological computations such as hydrodynamic 
models of the cloud. 
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SUPPLEMENTAL MATERIAL 
THOMAS-FERMI CONDITIONS 

We obtain the low-temperature EOS and central chem¬ 
ical potential /rg from a TF fit to the central region of 
the cloud. This method has the advantage of being in¬ 
dependent of a particular model for the EOS. It applies 
when the result is independent of temperature and fitting 
range, as we discuss in this section. 

We first motivate the linear scaling ep (/r) = c /I at zero 
temperature. In particular, we show that the density 
does not scale like /i, although ^ is positive on the Fermi 
side of the crossover. Due to the nonvanishing binding 
energy ep > 0 for all values of the interaction strength, 
the energy density (internal energy per volume) of the 
homogeneous system can be written as 

e(n) = —nd£B -I- 5c(n) = + <5e(n), (1) 

where = n/2 is the density of pairs. The chemical 
potential is found from 

de(n) £b , d(5e(n) , . . , 

The first term is independent of the density and thus 
does not contribute to the compressibility of the sample. 
This also uniquely defines the second term 5^{n). We 
write S^{n) = e-pjc, where the dimensionless constant c 
depends on the crossover parameter ln(fcFa 2 D)- Inverting 
the relation for ep we arrive at 

£p = y) = c il. (3) 

This is the claimed scaling and we further conclude that 
= /2(n). The quantity jKji) thus encodes the dif¬ 
ference between the interacting many-body system and 
a noninteracting gas of dimers. 

To justify the TF model at nonzero temperature we 
express the full EOS of the system in terms of the Fermi 
energy as 


£f(^,T, a2D) =/I • (4) 


with a dimensionless function g{x,y). For large /3/i 
(which corresponds to a large PSD) we can expand this 
expression according to 


£F(/r,r,a2D) = y(^go{l3eB) + ^gii/Sen) + ■ ■ - ^ (5) 


with go{y) = g{oo,y) and gi{y) = dxg{oo,y). The lead¬ 
ing term corresponds to the TF model with c = go(/3£B)- 
Below we demonstrate that the dependence of g{j3eB) 
on the binding energy is typically logarithmic in 2D. As a 
consequence, c is to a good approximation temperature- 
independent for a fixed magnetic field B. We can employ 
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FIG. 4: Temperature dependence of the TF parameter c on 
the Bose side (upper panel) and in the crossover regime (lower 
panel). We find c to be approximately temperatnre indepen¬ 
dent for an increase of T/Tp by 40% and 100% on the Bose 
and Fermi side of the crossover. 


this property to formulate a measure for the breakdown 
of the TF approximation; The validity of the expansion 
requires c to be independent of temperature. In Fig. 
0] we show the results of the fit of c in the Bose and 
crossover regime as a function of T/Tf . Each data point 
corresponds to the average of approximately 30 shots. 
We observe c to be approximately temperature indepen¬ 
dent for an increase of T/Tf by 40%(100%) on the Bose 
(Fermi) side. 

To derive a second condition for the applicability of 
the TF fit, we note that a strictly linear dependence 
£f oc /i would result in a constant c which is indepen¬ 
dent of the chosen fitting range. If the function is 

slightly nonlinear, the effective c will vary with a change 
in the fitting range. To quantify this influence we fit our 
data in three ranges of densities: (I) [0.4, 0.8]npeak, (H) 
[0.5, l]npeak, (III) [0.7, l]npeak, where npeak is the peak 
density. For a function £f(/I) which is slightly bending 
upwards for increasing /2, the fitted value for c would be 
smallest in (I) and largest in (III). We find the variation 
of c with the fitting range to be compatible within the 
errors in all three intervals. We observe the deviations 
to be smaller on the Bose side (5-10%) and to become 
larger as we approach the crossover regime and Fermi 
side (10-15%). There is another contribution to this ef¬ 
fect from the influence of the third dimension. As shown 
below, there is a //-dependence of £b, which results in 
ln(/3£B) = ln(/l£g^) -I- Aw{j^). Still, we find that both 
Cl and Cm agree with cu within the errors. 















To obtain a reliable estimate for the true value we de¬ 
fine the mean as 

( 6 ) 

where either C = /io,TF,C' = /r/ep, or C = ln(fcFa2D), 
and Cl/ll is the average of the observable over approxi¬ 
mately 30 shots with the TF model applied to the fitting 
range I/II. The corresponding statistical error is given 

by AC = yj\ {ACl + ACj^j) « ACi/n. The fitting range 
dependence of the TF model results in a systematic error 
in the determination of the low-temperature EOS. The 
further error due to systematic uncertainties is discussed 
below. 

Whereas we determine the low-temperature EOS from 
the measurements at the lowest attainable temperatures 
and estimate the error as outlined in the previous para¬ 
graph, it is important to estimate the temperature effect 
on the observable 1/c. Eor this we quantify the constant 
behavior of the function c{T/Tp) for small T/Tp (i.e. 
large phase space densities). We estimate the size of the 
plateau where c is constant by determining the plateau’s 
left and right boundaries (T/Tp)ynin < (T/TF)max from 
Fig. m 

For a fixed temperature T we define 1/c as in Eq. 

8 according to (1/c)t = to average the 

;ing range dependence and obtain the corresponding 

error A(1/c)t = . We then average 

over the data points for those temperatures within 
the plateau, ^ Ti < T^x- Let us assume that this 
comprises N data points. We then define the tempera¬ 
ture average as (l/c)temp = jf /^)Ti ■ The observ¬ 

able (1/c) temp has two main error sources: The fitting 
range dependence, and the spread of the values of (1/c)t 
with T. We generalize the fitting range error according to 

^(^/^)i‘ange — \/w and define the error 

due to averaging over temperatures as the standard devi- 

ation A(l/c)avg = jf Efci[(l/c)T. - (l/c)temp]^. The 
to tal error of (l/c)temp is then given by A(l/c)temp = 
yj A(l/c)?ange + A(l/c)2^g. We find that A(l/c)avg is 
significantly smaller than A(l/c)range and A(l/c)temp — 
A(l/c)T„i„. 

With the definitions of the previous paragraph we are 
able to compare the value of (l/c)T„i„ at the lowest 
attainable temperature Tmm and the average (1/c) temp 
over the plateau Tmin < T < T^ax- We show the 
results of this comparison in Tab. The ratio t = 
(T/TF)max/(T/rF)min — 1 gives the relative increase in 
T/Tp. The displayed values of 1/c thus provide a quan¬ 
titative statement on the size of the plateau. We find 
that the central value (l/c)temp is at most slightly shifted 
with respect to (l/c)T„i„ that A(l/c)avg is small such 
that A(l/c)temp — ^(l/c)T„i„- This indicates that the 
spread in temperatures of (1/c)t is small. Eurthermore, 


R[G] 

f _ ('^/^F)max 1 

(l/c)T„i„ 

(1/ c)temp 

692 

40% 

0.024(2) 

0.024(2) 

732 

49% 

0.046(4) 

0.045(5) 

782 

44% 

0.12(1) 

0.12(1) 

812 

38% 

0.21(1) 

0.22(2) 

832 

43% 

0.28(3) 

0.27(3) 

852 

52% 

0.38(4) 

0.38(4) 

892 

97% 

0.48(5) 

0.49(6) 

922 

93% 

0.54(6) 

0.53(6) 

952 

93% 

0.59(6) 

0.58(7) 

982 

96% 

0.59(6) 

0.59(7) 

1042 

106% 

0.67(8) 

0.64(7) 


TABLE 1: Thomas-Fermi fit. The slope c is approximately 
constant for small values of T/Tp in the interval (T/rF)min < 
T < (r/rF)max, see Fig. We show the size of this in¬ 
terval (plateau) in terms of the relative increase t in TjTp. 
We estimate the variation within this temperature-plateau by 
comparing the low-temperature EOS Ijc^ fi/ep at the low¬ 
est attainable temperature Tmin with the one averaged over 
the plateau. We label these and (l/c)temp, respec¬ 

tively. We observe the result of the temperature averaging 
within the plateau to hardly modify the value of 1/c and its 
error. The observable (l/c)T„i„ coincides with 1/c shown in 
Fig. 1 of the main text and displayed in Tab. |IV| below. 

the error due to the fitting range dependence is also only 
slightly changed due to averaging over the temperature- 
interval Tmin < T < Tmax- We conclude that the low- 
temperature EOS is faithfully captured by using either 
(l/c)T„i„ or (l/c)temp- For the remaining analysis we 
employ the former choice. 

TEMPERATURE DETERMINATION 

In order to determine the temperature of the sample we 
fit the outer region of the cloud to a reference EOS. Here 
we quantitatively compare different approaches and mo¬ 
tivate the particular choices employed in the main text. 

As our strategy for obtaining the EOS relies on a shot- 
by-shot determination of T and /tq, we cannot use the 
temperature determination from the tail of the momen¬ 
tum distribution from Ref. [7] , as it only gives the average 
over sets of shots. In fact, we find relatively large fluctu¬ 
ations from shot to shot for the temperature. Therefore, 
we display the median and the median deviation in the 
tables below. 

We first discuss the fitting ranges applied for the tem¬ 
perature fits. In each density profile we find that the 
local density fluctuations Aria ^ are below a thresh¬ 
old ~ 0.02^m“^ determined by the temperature of 
the gas as well as the imaging setup. The value of 
sets the density scale where the temperature fit can be 
applied. Since the density fluctuations average to zero. 
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{Aricr) = 0 , the density profile in the wings can safely be 
fitted also below the threshold, With tem¬ 

peratures T ~ 60 — 25nK on the Bose and Fermi sides we 
find ~ 0.2 — 0.4. 

To test for the fitting range dependence of the tem¬ 
perature fits of n(/i, T) we apply each reference EOS dis¬ 
cussed below to two fit regimes: These are (A) the den¬ 
sity range from O.lripeak to zero and (B) the /r-interval 
A/i = ±20nK centered at We further vary the size 

of /i-bins in the range 1 — 5nK. In all cases we obtain 
results for T and fio in close agreement with each other. 
For the following we thus restrict the discussion to fitting 
range (A) and a bin size of approximately InK. 

On the Bose side of the crossover, i.e., for magnetic 
fields B[G] = 692,732,782, the low density region for 
small g is given by the HF formula 

HF: ndAd = - ln(l - (7) 

For 692G we have g = 0.60 and HF theory should work 
reasonably well as has been demonstrated for bosons with 
g < 0.5 |42]. For larger couplings the HF approximation 
breaks down. A simpler formula for the Bose limit is pro¬ 
vided by a bosonic Boltzmann formula ridA^ = In 

the QMC calculations of the Bose gas we find the den¬ 
sity wings of the interacting Bose gas at 782G {g = 2.75) 
to be close to the bosonic Boltzmann formula even for 
moderate densities where deviations would be expected. 
The Boltzmann EOS is thus a trustworthy reference for 
the wings in this case. For the intermediate interaction 
strength g = 1.07 at a magnetic field of 732G, both HF 
and bosonic Boltzmann formula give the same tempera¬ 
ture, so we can use the Boltzmann result. 

A continuous interpolation between the bosonic and 
fermionic limits is provided by the second order virial 
expansion 

V: n,A^ = ln(I + e'^^) + 262e2'^'^, ( 8 ) 


In order to estimate how much the second order virial 
expansion deviates from the true EOS we define the ef¬ 
fective temperature exponent 


aesifi.T) = - 


dln(nCTA|.) 


dln^; 


1 d{ncrXx) 
n„\% d(/3^) 


( 10 ) 


For an EOS given by n^X'^ = z“ we have aeff = a for 
all /i and T. Eor instance, the fermion Boltzmann gas 
corresponds to cteff = 1 - In contrast, for large & 2 , the 
second order virial expansion gives = 2 , in accor¬ 
dance with a bosonic Boltzmann gas. In general, T) 
will not be constant. Here we are interested in its value 
for /3/i G I = [—2,0], which is the typical fitting range 
for the temperatures in our analysis. 

In Fig. we show a^e obtained from a numerical 
derivative of the LW and QMG data for the PSD. We 
observe that aeff deviates from the second order virial ex¬ 
pansion in the interval I for small values of , where 
it is close the fermionic Boltzmann limit a = 1. For 
/3eB ^ 5j the bosonic Boltzmann approximation a^s ~ 2 
is still a good account of the low-density region. The 
deviation from the second order virial expansion formula 
becomes relevant for /Seb ~ 2 , where aeff interpolates 
between the second order virial expansion value and 1 in 
the interval I. 

From this analysis we get a better understanding of the 
transition from a Bose to a Fermi gas. Indeed, when ap¬ 
plying a Boltzmann formula to the low-density region of 
the density distribution, we need to specify whether the 
degrees of freedom of the gas are bosonic or fermionic, 
i.e., whether the particles that behave classically are 
dimers or atoms. From the proximity of aeff to 2 for 
/3 £b ^ 5 we conclude that a good definition of the Boltz¬ 
mann formula for our parameters is given by 


B: Ua = 


{B < 832), 
(B > 832). 


( 11 ) 


where 


62 = 



exp{-eV(27r)} 

TT^ -I- (s — ln(27r/3EB))^ 


(9) 


is the (interaction induced correction to the) second virial 
coefficient (3^ [^. Equation Q is exact for small 
fermion fugacity z = in the whole 2D BEG-BGS 
crossover. Its applicability for moderate values of z, 
however, is not guaranteed. For large /Jeb, z ^ 0 and 
62 — such that Eq. ^ approaches the bosonic 

Boltzmann formula. As /3 eb becomes smaller, the sec¬ 
ond order virial expansion approaches the ideal Eermi gas 
formula. A comparison to LW and QMG theory, both of 
which should be reliable for sufficiently small /3 £b, indi¬ 
cates that the second order virial expansion is likely to 
overestimate the temperature of the sample, resulting in 
a curve n/no vs. 13fi which is too steep for small /3/r. 


The discontinuity in this definition leads to a jump in 
temperature obtained from the Boltzmann fit, see Tab. 
[TTl Note that we choose fi for the fermions instead of /i in 
Eq. ( [TT| ). Hence we have ncrA|. oc 5“ in the Boltzmann 
limit with z = This finding is motivated by the 

values obtained with the fit as summarized in Tab. [ml 
but does not affect our data analysis which is independent 
of the Boltzmann chemical potential. 

Note that a = 1 is o priori not a lower bound, and 
this is also not suggested by the theoretical data for 
even smaller /3 eb- Still, for the interaction parameters 
5 ^ /3 £b ^ 0.5 relevant to Fig. 3 of the main text, the 
best available theoretical predictions indicate that aeff 
interpolates between the second order virial result and 
the fermion Boltzmann formula a = 1 in the interval I. 

Denoting the temperature obtained from a fit of Eq. 
([^ to the low-density region by Ty and the one using Eq. 
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FIG. 5: The effective temperature exponent aefi defined in 
Eq. (101 indicates the deviation of the system from a bosonic 
and fermionic Boltzmann gas (aefi = 2 and Qeff = 1), shown 
as dotdashed lines in both panels, for /3/i € X = [—2, 0]. Fur¬ 
thermore, Qeff has been obtained from a numerical derivative 
of the data for the PSD from LW theory ([23]) and QMC 
( [24] 1. This is not meant to be a precise characterization of 
the corresponding EOS, but should rather indicate the over¬ 
all trend. Upper panel. The upper panel shows the predic¬ 
tions for Peb. = 5 from LW theory (green solid line) and for 
Peb = 3 from QMC (red dashed line). The dotted lines are 
the corresponding curves from second order virial expansion. 
We observe q = 2 to be a good approximation of the up¬ 
per curve in T. Lower panel. We compare LW (solid lines), 
QMC (dashed lines) and the virial expansion (dotted line) 
for Peb = 2 (upper curves) and Peb = 0.5 (lower curves). 
The data suggest to interpolate between the virial curves and 
a = 1 in the fitting range X. 


0 by Tb, we thus conclude that 

y _ jy + 


( 12 ) 


is a good estimate for the temperature of the sample for 
Peb ^ 0.5: On the Bose side, Ty — Tb and the average is 
trivial. On the Fermi side, our analysis of Oeff indicates 
that Ty and Tb are upper and lower bounds to the ac¬ 
tual temperature. We apply the temperature definition 


in Eq. (12) for the magnetic fields 732 < i3[G] < 1042. 
The systematic error for B[G] > 892 where aes becomes 
smaller than 1 is not relevant for the results on the EOS 


B[G] 

Ttof 

Thf 

Tb 

Tv 

T 

692 

55(18) 

61(8) 

a = 2 

51(4) 

51(4) 

61(8) 

732 

55(20) 

52(5) 

a = 2 

49(5) 

49(5) 

49(5) 

782 

46(14) 

36(3) 

a = 2 

47(5) 

47(5) 

47(5) 

812 

57(9) 

- 

Q = 2 

45(7) 

45(7) 

45(7) 

832 

59(13) 

- 

a = 2 

40(8) 

37(6) 

39(7) 

852 

58(14) 

- 

a = 1 

22(4) 

36(5) 

29(4) 

892 

71(19) 

- 

a = 1 

23(3) 

33(3) 

28(3) 

922 

67(19) 

- 

a = 1 

23(4) 

30(5) 

27(5) 

952 

57(11) 

- 

a = 1 

20(3) 

26(4) 

23(3) 

982 

67(26) 

- 

a = 1 

23(5) 

28(5) 

26(5) 

1042 

45(15) 

- 

a = 1 

23(6) 

26(6) 

25(6) 


TABLE IL Lowest temperatures (in nK) attained in our ex¬ 
periment. We compare the results from fitting with different 
reference EOSs. The first column corresponds to the temper¬ 
ature obtained from a Boltzmann fit of the tail of the mo¬ 
mentum distribution after time-of-flight (TOF) [7]. The cor¬ 
responding values are the mean and standard deviation after 
averaging over ~ 30 shots. The remaining columns show the 
median and median deviation of fitting the low density region 
of the in-situ density profile to the Hartree-Fock temperature 
Thf 0 , Boltzmann temperature Tb with a (111, and second 
order virial expansion Ty ([^. The last column displays our 
best estimate for T as discussed in the text. Note that we de¬ 
termine the temperature for each shot individually, such that 
the mean values and their error displayed here only serve as 
an orientation. 


B[G] 

Po.TF 

Ao.hf 

h0,B 

Uo,v + ^ 

ho,v 

692 

34(3) 

39(3) 

49(5) 

- 

- 

732 

51(4) 

60(5) 

65(5) 

- 

- 

782 

78(6) 

91(6) 

75(8) 

- 

- 

812 

110(7) 

- 

100(6) 

110(10) 

-180(10) 

832 

130(7) 

- 

120(10) 

120(6) 

40(6) 

852 

160(8) 

- 

170(6) 

140(8) 

110(8) 

892 

200(7) 

- 

200(7) 

180(9) 

170(9) 

922 

210(10) 

- 

200(10) 

190(10) 

190(10) 

952 

190(10) 

- 

180(10) 

170(10) 

170(10) 

982 

210(10) 

- 

200(7) 

190(6) 

190(6) 

1042 

200(10) 

- 

180(20) 

180(10) 

180(10) 


TABLE III: Central chemical potentials in nK. The first col¬ 
umn shows the mean and standard deviation for the Thomas- 
Fermi (TF) value Po.tf which is used in the main text. From 
the HF-, B-, and V-fits of Tab. |II|we obtain effective chemical 
potentials, which are close to the TF result. We show their 
median and median deviation values. The 2D binding energy 
eb which relates /io and po is given in Tab. |IV| below. 


presented in our article. Only for B = 692G we apply 
the better estimate from HE theory as a Boltzmann fit 
is clearly disfavoured by the tail of the density profile. 

In Tab. [Tljwe summarize the results of the temperature 
fit using the reference formulas presented in this section. 
We indeed find the bosonic Boltzmann values to be close 























































11 


to the result from the virial expansion for B < 832G. For 
larger magnetic fields we confirm the expectation that 
the virial fit is consistently above the fermionic Boltz¬ 
mann result, although both approach each other. For 
B > 892G the deviation between Ty and T-q is lOnK or 
less. We note that Ty gives a smooth temperature varia¬ 
tion as we cross the Feshbach resonance from low to high 
magnetic fields. In contrast, Tb jumps due to the discon¬ 
tinuous definition of a in Eq. We conclude that the 

smooth phenomenological interpolation for a employed 
in [7] is justified by the second order virial expansion. 

Further note that the decrease in temperature from the 
Bose to the Fermi side of the crossover is approximately 
by a factor of two. This can be explained by assuming 
an adiabatic ramp across the Feshbach resonance when 
changing the interaction strength. Since the degrees of 
freedom double at constant entropy, the temperature is 
reduced by a factor of two. We find a good agreement 
with the temperature obtained from the large momentum 
tail of the momentum distribution after TOF on the Bose 
side, but find substantial differences on the Fermi side, 
as has already been discussed in Ref. [7] . For the present 
analysis we obtain the temperature from fits of the in-situ 
density profiles as this gives an overall consistent picture. 

From the fit of the outer regions of the cloud we ob¬ 
tain, as a by-product, an estimate for the central chem¬ 
ical potential /tq. We summarize the corresponding val¬ 
ues in Tab. Ill We also show the chemical potential 
determined by the TF fit of the high-density central re¬ 
gion. We find good agreement within the errors from 
the different methods. This indicates that the reference 
EOSs introduced in this section are well-applicable to the 
wings. Furthermore, deviations in fiQ obtained from ap¬ 
proximate reference EOSs does not affect our data anal¬ 
ysis based on /io.TF, which we have shown to be robust. 


2D SCATTERING PHYSICS 


corrections are important if /(fco) becomes of order unity, 
where ko is a typical momentum of scattering particles. 

The scattering amplitude for the 3D system which is 
strongly harmonically confined in the z-direction with an 
oscillator length is given by mm 


/q2D(fc) 


Att 


where the function w{^) reads 





J 


E 


(2j-l)!! 

(2j)!! 


ln(j - C - iO) 


(14) 


(15) 


For low density ^ <C 1 we have w{^) ~ rciim(C) = 
— A) + iTT with A = 0.905. We define the filling 

correction 


Au;(^) = w{C) - wiim(C) > 0. (16) 


For ^ < 0 we define Aw(^) = 0. If is sufficiently 
small (tight confinement), we have Aw « 0, and fq 2 B{k) 
approaches the form of f 2 T>(k) with a 2 D given by 



(17) 


Neglecting Aw is justified if the relevant momentum 
scale of scattering, kq, satisfies fco^z ^1. In a many- 
body setting, the corresponding momentum scale is given 
by hko = y/2Mfl [6], thus kQ£ll2 = jl/huz- If is not 
small, the effective 2D scattering length of the many- 
body system receives a fcof^-dependent correction. The 
latter is obtained from keeping the higher order terms in f 
in the denominator of Eq. (14) when equating f 2 B{ko) = 
/q 2 D(Ao)- We then arrive at the more general formula 


Due to the tight confinement along the z-direction, the 
scattering physics of the 3D atomic gas can be mapped to 
those of a 2D gas. Here we summarize the formulas which 
are required to translate the 3D scattering length a 3 u{B) 
across the Feshbach resonance at 832.2G m to the 2D 
scattering length characterizing the quasi-2D system. For 
a more detailed discussion of similar aspects see also Refs. 

01711131 . 

Scattering in 2D is peculiar in the sense that the scale 
dependence of interaction corrections is typically loga¬ 
rithmic. This is reflected in the universal low-momentum 
behavior of the 2D scattering amplitude given by 


a2D 



\l\w{klAj2) 


(18) 


For the lowest temperatures attained in the experiment, ^ 
is irrelevant in the bosonic limit, but reaches ^ ~ 0.7—0.9 
for the highest magnetic fields on the Fermi side. The 
shift of ln(/cFa 2 D) with respect to ln(fcFa 2 B), however, is 
only — ^Aw = —(0.4 — 0.5) in this case. 

We define the 2D binding energy as 


- 


(19) 


/2d(A:) 


47r 

— \n{k?a\j^) + iTT 


(13) 


for fc —>■ 0, where k is the relative momentum of two col¬ 
liding particles in the center-of-mass frame. Interaction 


Analogously, with /{M[a^^]'^) we have £b = 

The effect of Aw appears now exponentiated. 
However, the dependence of the PSD on /Ieb is only log¬ 
arithmic. Note that eb does in general not coincide with 
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the energy Eg of the confinement induced bound state [T] 
given by the solution of IzIa-iVi = /i(£b/(^z)) 


potential V (r). For the purpose of the present analysis 
we assume 


Jo ^ 




^(1 _ e-2“)/(2'u) 
( 0 ) 


( 20 ) 


We compare the values of Eb , £b > ^b magnetic 
fields 692 - 1042G in Tab. ^ 


SYSTEMATIC UNCERTAINTIES 

In the main text we display the statistical errors of our 
determination of the EOS. Due to systematic uncertain¬ 
ties in the measurement of the density profiles n(r), from 
which the EOS is constructed with the EDA, we have fur¬ 
ther systematic errors which are summarized here. Note 
that these systematic errors do not cover the deviations 
which arise from our choices of the fit function to de¬ 
termine the thermodynamic parameters T and /tq- The 
latter can hardly be quantified at this point. The system¬ 
atic errors of the experiment have already been discussed 
in detail in Ref. [7] . 

Absorption imaging. The planar density profile 
n(r) is constructed from the (column integrated) opti¬ 
cal density OD(r) after absorption imaging along the z- 
direction. At zero detuning we have 


n(r)(TQ = OD(r) -|- 


7o(r0 


(1 


_ „-OD(f 


'), ( 21 ) 


where Ug is the effective absorption cross section, Jg the 
initial intensity of the probe beam before the atomic 
cloud, and is the effective saturation density. The 
experiments are performed for lo/Igat = 0-97^°;Jg. This 
leads to a systematic error in the peak density npeak^ 4 %- 
Furthermore, due to the nonzero binding energy, the op¬ 
tical transition frequencies of dimers are shifted, which 
leads to a reduced absorption cross section CTq on the 
Bose side. We compensate for this effect for i?[G] = 
692,732,782 by means of a rescaling in CTq. The sys¬ 
tematic uncertainty in this rescaling factor leads to an 
additional systematic error in the peak density ripeak of 
at most 8%. 

Atoms in noncentral pancakes. Our experiment 
realizes the planar gas in a pancake shaped geometry 
at z = 0. Therein the pancake is realized as a standing- 
wave optical dipole trap from interfering two laser beams 
at a small angle each. In this way, further pancakes at 
z ^ 0 are created at a distance ~ 4^m from the central 
pancake. A nonvanishing population of the noncentral 
pancakes leads to an overestimation of the planar density. 
Assuming the gas in the higher pancakes to be a thermal 
gas we find that the density is overestimated at most by 
10 %. 

Trap parameters. The application of the EDA re¬ 
quires an accurate knowledge of the in-plane trapping 


V{f^ = ^{colx'^ + + caxx'^ + C 4yy^ + c^xyX^y 


^y\ 

( 22 ) 


where the potential is centered at ro = 0 , and oJx = ujy 
[7] .The harmonic trapping frequencies uJx/y are known to 
an accuracy of 0.4%, hence their uncertainty can safely 
be neglected in the following. The quartic terms C 4 are 
only known to a precision of 25%, but they are small 
and hardly contribute to the density region of interested. 
Higher order terms in the expansion are not relevant as 
the gas does not extend to regions so far from the center. 

Magnification. The magnification of the imaging 
system relates the pixels of the imaged density profile 
to in n^f). The magnification factor is known to 

an uncertainty of 3%, which results in an error of 6 % in 
V (r), where the distance appears squared. 

Binding energy. The determination of the binding 
energy from eb = requires knowledge of 

030 (.B) and G to compute Eg^ Both quantities, however, 
are known very well in our experiment and their error can 
be neglected [7]. On the other hand, for large /Ig, the de¬ 
viation due to the exponential of Aw can be substantial. 
We have hwz/k^ = 265nK. Assuming a deviation of )io 
by ±10% from the mean values /2o,tf displayed in Tab. 
Ill we find the systematic error in eb to be ±4%. 

The listed systematic uncertainties affect the thermo¬ 
dynamic variables n, y,, and T. The density is influenced 
by the uncertainties in the absorption imaging. This 
leads to an uncertainty for B < 782G (Bose side) 

and for B > 812G (crossover regime and Fermi 

side). The chemical potential is derived from the EDA as¬ 
sumption /i(f^ = yo—V (r) and hence inherits the system¬ 
atic uncertainties in U(F). Those result from the uncer¬ 
tainty in the trapping frequencies and the magnification 
of the imaging system, giving . The temperature is 
fitted in the low-density region with a reference EOS of 
the form = X^‘^v[j3yQ — /3U(r)), where v{x) is a di¬ 
mensionless function of a; = I3y. The uncertainties of the 
density influence both and yo- For instance, in the 
Boltzmann case we have Uo- = _ (j^-PVir)^ 

such that density uncertainties are absorbed in the ir¬ 
relevant prefactor C. We conclude that the systematic 
uncertainty in T is the same as the one for y, i.e. ± 6 %. 


MAXIMUM OF THE EQUATION OF STATE 

The EOS at nonzero temperature expressed by the 
function n/ng = h{f3y, Psb) exhibits a maximum as a 
function of Py for fixed /3 eb. Here we quantify the loca¬ 
tion and height of this maximum from our experimental 
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B [G] 


020 [fim] 

4°^[nK] 

Eb [nK] 

4[nK] 

Ao[nK] 

Ep [nK] 

ln(A:Fa2D) 

Mo/EF 

1/c 

692 

0.000137 

0.000128 

4.31 xlO® 

4.96 xlO® 

13600. 

33(4) 

1180(123) 

-7.27(5) (/«) 

0.028(4) (+^) 

0.024(2) (+^) 

732 

0.00703 

0.00624 

1.63 xlO® 

2.07 xlO® 

4330. 

50(6) 

951(132) 

-3.49(7) (/«) 

0.052(10) (/«) 

0.046(4) (/^) 

782 

0.147 

0.12 

3730. 

5610. 

766. 

78(7) 

590(44) 

-0.76(3) (/«) 

0.13(2) (/"J 

0.12(1)(/2) 

812 

0.517 

0.376 

302. 

569. 

191. 

113(9) 

469(36) 

0.29(3) (+®) 

0.24(3) 

0.21(2)(/3) 

832 

1.02 

0.699 

77.3 

165. 

65.5 

130(10) 

417(45) 

0.86(5) (+®) 

0.31(4) (/^) 

0.28(3) (/^) 

852 

1.83 

1.15 

24. 

61.5 

22.6 

154(13) 

371(41) 

1.31(5)(+«) 

0.41(6) («) 

0.38(4) (/5) 

892 

4.73 

2.45 

3.61 

13.5 

3.58 

193(16) 

366(35) 

2.10(8) (+®) 

0.53(7) (/^) 

0.48(5) (/^) 

922 

8.29 

4.01 

1.17 

5.00 

1.17 

204(18) 

341(40) 

2.58(8) (+®) 

0.60(9) (/*) 

0.54(6) (/^) 

952 

13.2 

6.96 

0.46 

1.66 

0.459 

191(19) 

295(37) 

3.02(7) (+®) 

0.65(10) (+®) 

0.59(6) (/«) 

982 

19.6 

9.48 

0.209 

0.897 

0.209 

205(14) 

319(25) 

3.40(8) (+®) 

0.64(7) (/®) 

0.59(6) (/«) 

1042 

36.5 

18.4 

0.0605 

0.237 

0.0605 

197(21) 

268(33) 

3.96(10) (/®) 

0.74(12) (+^0) 

0.67(8) (/®) 


TABLE IV: Scattering parameters and low temperature EOS as shown in Fig. 1. We show the result of averaging the 
observables over approximately 30 shots for each magnetic held. The error is shown as (stat.) , where the statistical error 

is given by the standard deviation and the systematic error results from the systematic uncertainties discussed in this section. 
Note that the 2D binding energy eb is generally larger than the quasi-2D universal dimer energy Eb- 
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FIG. 6: Comparison of the measured location and height of 
the maximum of the curve n/no (red, “Expt”) with theo¬ 
retical predictions from LW ([23], green) and QMC ( [24] . 
orange) theory. In the upper panel we show the location 
of the maximum (/9^)max for hxed /3eb. In addition to 
the expected linear scaling with —/1 eb/ 2 we find a sublead¬ 
ing correction of « ln(2). The solid blue line is the curve 
(/3^)niax = —/3 eb/2 + ln(2). In the lower panel we show the 
height (n/no)max of the maximum as a function of /3eb. The 
blue solid curves is the reference (n/no)max = 


data and compare to theoretical predictions from QMC 
and LW calculations. 

Our data is consistent with a maximum of height 
(n/no)max ^ 26^®=^/^ at (/3Ai)max ^ -^+ln(2). Writing 
X = (3fj, and y = (n/no), the maxima thus all approxi¬ 
mately lie on the curve y = 4e““. We demonstrate this 
behavior in Fig. where we also compare to theory. 





































